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۲ ,ر)مططعه]۷ ب‌مططمها له تعمهجنصتا تعهموظ رععصمنگ آم‌نامهمطها له مج 1 
میم ]تمصصم )12۲26 2۵ 


7۲صقحظ. نود 0معصمتتص مق اقطا قصمحصمصمطم عمجم مه فعفهمودنا. ,ءاعهافواظ 
8 3۵6۵6 مصتولتاای جمل امه احمهعوص بچ مدز متام مهد لهمزامصمطمما . ,قمامه] 
وعوهعنل 0صهماعملصه میا ومرامط مصتامل‌ما .قعمم‌معور مققععتل فطع م فمععله تتمطاه 220 
۶ واممج اهعتامصمطها محصقطه امجمم مه فوع عتاععه متمصظد 0فصق 0صنق اتمه 
هه مه مطاتععمل اقطه مصمتمنوه همم حدم 0عصقط «الفنامده مته فععفعوتل 
0 صق فلمل‌مصط معمط1 .عصصا من مهوت مطل طنه ماممم گم صمتعلدومص مط ۵ 
احمصاهما ۵۶ فعمصم۲تاعه مطا ماملعی رمعمع‌وتل فطع ۵۶ مقعنامی فطل 606 مت 121660 
۶ مصو) ,عصهععتل مطه صم مه مامتنعب و له مطه ماهعتافمط صنه رقفل0ظ۵ظ 
مقط (06مصظ فنط 1 ,مهد عو مط قوذ فع‌عهمعنل ۶ علمل‌مصط آمتامصصمطممه اعفمامرحصته فطل 
مط اعتعجر ما 0عفیها وا مضه فععهمعت میامتامعص گم تجناای مط صا مصمتام‌تاموه ۳100 
مط ماممت)امم۲صا مصع رصمتامصی ۲۵ ۵۶ وعمصمتا0ع له مطا ماععیي رم‌تصهلتوو ۵۶ معصامع 
6۰ ۱( ۵۴ موعیامم فطل دم متحقطه‌ه صمصصیبیط قق طمیاو فجمشه منامتنوه ۵۶ عم 
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مقاله علمی-ترویجی 


به سوی علوم ریاضی 


۵ ۴ص عمط 05۰/۱ 


مدلسازی ریاضی و تاثیر آن در شناخت و کنترل بیماری‌ها 


۲ ۱۶۰ 
مرجب» زاج 
دانشکده علوم ریاضی. دانشگاه فردوسی مشهد ایران 


جرمی.]تمصصم )م2 ۲ق1. 2۵ 


چکیده. بیماری‌ها پدیده های پیچیده ای هستند که تحت تأثیر غوامل متعددی قرار دارند. مدلسازی ریاضی 


ابزاری قدرتمند برای مطالعه این عوامل و تأثیر آنها بر روند بیماری است. در واقع مدلسازی به ما کمک کند 
تا بیماری‌ها را بهتر بفهمیم و راه های موثرتری برای کنترل آنها پیدا کنیم. مدل‌های ریاضی بیماری‌ها معمولاً 
بر اساس معادلات ریاضی هستند که رفتار جمعیت افراد مبتلا به بیماری را در طول زمان توصیف می‌کنند. 
این مدل‌ها می‌توانند برای پیش بینی روند بیماری» ارزیابی اثربخشی روش‌های درمانی و بررسی تأثیر عوامل 
مختلف بر بیماری استفاده شوند. یکی از ساده ترین مدل‌های ریاضی بیماری‌ها. مدل ۹115 است. این مدل 
کاربردهای گسترده ای در مطالعه بیماری‌های عفونی دارد و برای پیش بینی روند اپیدمی‌ها؛ ارزیابی اثربخشی 
واکسیناسیون و بررسی تأثیر عوامل مختلف مانند رفتار انسان بر روند بیماری استفاده می‌شود. 


۱ مقدمه 


بیماری‌ها پدیده های پیچیده ای هستند که متاثر از عوامل متعددی پیش می آیند. یکی از ابزارهای قدرتمندی 
که می‌توان برای مطالعه و شناخت بیماری‌ها و عوامل رشد و شیوع انا سکار تفرگ مسارم وتا طی استت: 
مدل‌های ریاضی به‌طور گسترده در طول همه‌گیری‌های گذشته و حال مورد استفاده قرارگرفته‌اند تا نگاهی دقیق‌تر 
به پویایی انتشار بیماری داشته باشند و پیش‌بینی کنند که چه زمانی و به چه قیمتی بیماری مهار می‌شود. در 


0 ,34115 .جملامصتعععان بممزطاند متامصجمط) 2 2020 
واژگان کلیدی. مدلسازی ریاضیی بیماری عفونی» مدل چند محفظه ای» اپیدمی. 
تاریخ: دریافت ۱۴۰۲/۷/۱ بازنگری ۱۴۰۲/۷/۲ پذیرش ۱۴۰۲/۷/۲ انتشار برخط ۱۴۰۲/۷/۵ 
نحوه ارجا ع به اين مقاله: م.زاج» مدلسازی ریاضی و تاثیر آن در شناخت و کنترل بیماری‌ها . به سوی علوم ریاضی» ۳ (۰)۱۴۰۲ 
شماره ۰۱ ۰۷۲-۶۲ 


9دانشگاه فردوسی مشهد. 
۶۳ 


۶۴ م. زاج 
پاب اهمیت مدلسازی جا دارد متذکر شویم که این علم یکی از ابزارهای مهم در مبارژه با بیماری گوید ۱۹ بوده 
و هست و به ما کمک می‌کند تا با تجمیع تحقیقات موجود در این زمینه» چالش‌ها و شکاف‌های تحقیقاتی را 
هر چه بهتر شناسایی کنیم و بتوانيماقدامات لازم را بای کنترل اپیدمی انجام دهیم. مدلسازی ریاضی معمول 
با استفاده از معادلات ریاضی برای توصیف تعاملات بین عوامل مختلف و تأثیرگذار بر یک بیماری ایجاد 
می‌شوند که اين معادلات را عموما به عنوان یک سیستم دینامیکی در نظر می‌گیریم. در واقع» سیستم‌های 
دینامیکی مجموعه ای از معادلات ریاضی هستند که رفتار یک یشیم را در طول زمان توصیف می‌کنند. این 
سیستم‌های دینامیکی می‌توانند برای پیش بینی روند بیماری» تأثیر عوامل مختلف بر بیماری و بررسی راه های 
مختلف درمانی استفاده شوند. مدل‌های ریاضی می‌توانند برای پیش بینی روند بیماری در اینده استفاده شوند. 
این تسشن منگ‌ها مق اک رای برنا مه ریوی ,هر یه بیقر برای مراقیتهای بهداشقی و افتامات کل تایه 
شوند. همچنین. می‌توانند برای بررسی تأثیر عوامل مختلف بر بیماری استفاده شوند. این عوامل می‌توانند 
شامل عوامل ژنتیکی» محیطی و رفتاری باشند. سرانجام با توجه به تمام اطلاعاتی که اين سیستم‌های دینامیکی 
و اتشتارمان فرار مه دنه تزا برای ارفا ای شلات رای اسظانه قوف این تفای 
مي‌توانتد. شامل واکستاسیزن فرمان دارویی و اقدامأث کشرلی باشند: 

یکی از این مدل‌های ریاضی. مدل‌های چند محفظه ای هستند که برای توصیف رفتارهای بیماری‌های 
واگیردار و عفونی مانند آبله ایدن مالاری کرونا و ... بکار می‌روند [۱۷]. مدل‌های چند محفظه‌ای از اوایل 
قرن بیستم نقش اصلی را در مدل‌سازی بیماری‌های عفونی ایفا کرده‌اند. این مدل‌ها به طور گسترده برای توصیف 
فرایند یک اپیدمی در حال گسترش استفاده می‌شوند. معمولا این مدل برای اپیدمی‌های کوتاه مدت بکار برده 
می‌شود. اپیدمی‌های کوتاه مدت اپیدمی‌هایی هستند که در مقایسه با طول عم مدت زمان آنها نسبتاً کوتاه 
است. این اپیدمی‌ها معمولاً به شکل شیوع ناگهانی بیماری هستند که بخش قابل توجهی از جمعیت یک منطقه 
را مبتلا می‌کنند و احتمالاً می‌کشند. به طور کلی در اپیدمی‌های کوتاه مدت. نرخ انتقال معمولاً زیاد است. این 
بدان معناست که افراد مستعد به سرعت به افراد آلوده تبدیل می‌شوند. همچنین؛ نرخ بهبودی یا مرگ و میر 
معمولاً زیاد است. این بدان معناست که افراد آلوده به سرعت بهبود می‌پابند يا می‌میرند. مدل‌های چند محفظه 
ای مانند مدل ۹1 علاوه بر پیش بینی روند اپیدمی. نشان دهد که اپیدمی چگونه گسترش می‌یابد و چگونه در 
نهایت فروکش می‌کند. 


۲ ساخت مدل بیماری عفونی 


بررسی قرار گرفتند. این مدل یک مدل ساده از انتشار بیماری است که جمعیت مفروض را به سه گروه مجزا 


۵۲ ویمتاهع کم ماطنامعععن6 


مدلسازی ریاضی و تاثیر آن در شناخت و کنترل بیماری‌ها ۶۵ 
تقسیم می‌کند. جمعیت افراد سالم و مستعد ابتلا به بیماری که می‌توانند از طریق تماس با افراد آلوده به بیماری 
مبتلا شوند را با 5 " نشان می دهند. جمعیت افراد آلوده به بیماری را با آ ۲ نشان می‌دهند که این افراد می‌توانند 
بیماری را به دیگران منتقل کنند و جمعیت افراد بهبود یافته را با 1 * نشان می‌دهند. مدل ٩11‏ با معادلات 


زیر توصیف می‌شود: 
(۱.۲( ]۰ - 68 < سب 


در این معادلات 8 نشان دهنده نرج انتقال بیماری و + نشان دهنده نرج بهبودی از بیماری فتاه در وأقع این 
یک مدل سه بخشی برای مطالعه چگونگی تکامل بیماری‌های عفونی در طول زمان در سطح جمعیت است. 


6 « 55 « سس 


شکل ۱: سیستم دینامیکی مدل سه محفظه ای مستعد-آلوده-بهبود یافته 


به طور خلاصه. این مدل وضعیت عفونت در جامعه را با (ّ) انتقال افراد مستعد به محفظه عفونی از طریق 
فرآیند انتقال و () انتقال افراد عفونی به محفظه 15 (چه مرده و چه بهبود یافته) از طریق فرآیند حذف توصیف 
می‌کند. همانطور که قبلا اشاره کردی» در یک زمان معین؛ کل جمعیت ۱ تحت یک مطالعه به سه بخش تقسیم 
می‌شود که با 5 ۰ ] و 1 نشان داده می‌شوند. و اندازه‌های آن‌ها در تساوی ۷ <- 0 + + 5 صدق می‌کند. 

اغلب اوقات. علت اصلی علاقه به چنین سیستمی, در مقادیر تابع در طول زمان نهفته است. با اين حال؛ 
ممکن است راه حل تحلیلی مناسبی برای چنین توابعی وجود نداشته باشد. به عنوان مثال, برای پاسخ به این 
سوال که تا پایان سال ۲۰۲۰ (یا هر زمان اینده) چند نفر به کووید-۱۹ مبتلا خواهند شد. نیاز به ماشین حسابی 
داریم که بتواند تعداد تجمعی موارد مستعد آلوده و حذف شده و بهبود یافته را از گذشته تا آینده محاسبه کند. 
متأسفانه, در واقعیت. توابع مرتبط با این ماشین حساب معمولاً غیرخطی هستند و فرم دقیق آنها را به سختی 
می‌توان به طور مستقیم مشخص کرد. در مقابل» مجموعه‌ای از معادلات دیفرانسیل معمولی (01917) به درک 

"اطناوعع:۹ 


10160۱0" 


13010 


۶۶ م زاج 


بهتر دینامیک انتقال بیماری (یعنی ویژگی‌های بیماری‌های عفونی) کمک می‌کند و ویژگی‌های کلیدی آنها را 
به‌طور راحت‌تری پدیدار می‌کند. در چنین حالتی» ری بیماری با یک مدل دستگاه دینامیکی معمولی ممکن 
است مطابقت داشته باشد زیرا معادله دیفرانسیل معمولی» یک معادله ریاضی است که رابطه بین یک یا چند 
متغیر و نرج تغییر آنها را توصیف می‌کند. این معادلات به طور گسترده در فیزیک» مهندسی و زیست شناسی 
برای مدلسازی سیستم‌های پیچیده استفاده می‌شوند. 

برای حل عددی معادلات دیفرانسیل معمولی با شرایط مرزی معین» می‌توان از روش‌های عددی مانند روش 
گسسته سازی اویلر یا روش تقریب رانگ-کوتا استفاده کرد [۲] و [۱۵] ۰ زیرا حل معادلات دیفرانسیل با 
تایط مری م ‏ ز اش عیرس ای کاس تفع پاش موش فا عفزی 
می‌توانند برای حل این معادلات به طور تقریبی استفاده شوند. 


مثال ۰۱.۲ مدل ۹15 را براي یک جمعیت فرضی با تعداد ثابت 100 -- 7۷ در نظر بگیرید» که از این تعداد 
٩‏ نفر مستعد. ۱ نفر آلوده و هیچ بهبود يافته یا حذف شده ای نداریم. انتقال بین محفظه هاء که در دستگاه 
معادلات ۱.۲ نوشته شده است» نشان دهنده حرکت جمعیت از یک محفظه به محفظه دیگر است. اگر مقادیر 
پارامترها را 0.5 < 8 و 0.2 < چ در نظر بگيریم, با استفاده از یکی از روش های عددی که قبلا اشاره کردیم» 
می‌توان رفتار جمعیت های محفظه های مختلف را طبق شکل ۲ با نمودار به خوبی توصیف کرد. 
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شکل ۲: جواب معادله دیفرانسیل معمولی مدل سه محفظه ای با استفاده از روش تقریب مرتبه چهارم رانگ کوتا 


۲ عدد بازتولید. بر اساس دو پارامتر 0 و7 دز یکت مدل ٩1۳‏ 6 نسبت # - 1۵ به عنوان عدد 
بازتولید پایه ای نامیده می‌شود که در واقع نشان دهنده تعداد افرادی است که به طور مستقیم ویروس را از یک 
فرد آلوده دریافت می‌کنند. 770 یک پارامتر مهم برای ارزیابی چگونگی گسترش یک بیماری در جمعیت است و 
به عوامل مختلفی مانند نرخ انتقال ویروس مدت زمان آلوده بودن افراد و جمعیت مستعد بستگی دارد. برای 
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مثال» 710 برای ویروس سرخک حدود ۱۸-۱۲ است. به این معنی که یک فرد مبتلا به سرخک به طور متوسط 
۲ تا ۱۸ نفر دیگر را آلوده می‌کند. در مقابل» 77 برای ویروس ۲1۲۷ حدود ۴ است. به اين معنی که یک 
فرد مبتلا به 117۷ به طور متوسط ۴ نفر دیگر را آلوده می‌کند. بنابراین» عدد بازتولید یک ویژگی کلیدی برای 
توصیف و مقایسه بیماری‌های عفونی ارئه می‌دهد. مراجع [۴] ۰ [۵]» [۸] و [۱۰] را ملاحظه کنید. با توجه 
به این اطلاعات می‌توان ادعا کرد که اگر 770 در یک بیماری بیشتر از ۱ باشد. انتظار می‌رود که آن بیماری به 
صورت اپیدمی گسترش یابد. در مقابل, اگر م71 یک بیماری کمتر از ۱ باشد انتظار می‌رود که آن بیماری از 
بین برود. در واقع بیشتر از ۱ بودن عدد م73 به این معنی است که هر فرد آلوده به طور متوسط بیش از یک نفر 
دیگر را الوده می‌کند. در این صورت. بیماری می‌تواند به سرعت در جمعیت گسترش یابد و یک اپیدمی ایجاد 
کند. در مقابل اگر 77 در یک بیماری کمتر از ۱ باشد. به این معنی است که هر فرد آلوده به طور متوسط کمتر 
از یک نفر دیگر را آلوده می‌کند. در این صورت. بیماری به کندی در جمعیت گسترش می یابد و در نهایت از 


بین می‌رود. 


۳ فرض‌ها و محدودیت‌ها در مدل 91۳۴ 


مانند هر مدل ریاضی دیگری» مدل ۹1 نیز دارای برخی فرض‌ها و محدودیت‌ها مانند شرایط مرزی است 
که باید برآورده کند. این محدودیت‌ها شرایطی را تعریف می‌کنند که در آن مدل 818 ممکن است در عمل قابل 
استفاده باشد. تمام فرض‌های کلیدی عبارتند و 


جمعیت بسته است. به اين معنی که هیچ فرد جدیدی به جمعیت اضافه نمی‌شود و هیچ فردی از 
جمعیث خارج نمی‌شود. مثلا ما در اینجا جمعیت را برابر با مقدار ثابت آ در نظر گرفته آیم. در 


و اه 

افراد در جمعیت به طور تصادفی با یکدیگر ملاقات می‌کنند. زیرا احتمال و درجه تعامل با یکدیگر در 
طول زمان ثابت می‌ماند. این یک فرض قوی از همگنی برای سیستم دینامیکی ۹11 است که توسط 
همان پارامترهای 8 و 7 کنترل می‌شود. فرض همگنی باعث می‌شود که مدل 5115 قابل محاسبه و 
تحلیل باشد. 

9 در یک فرد مستعد فقط از طریق عفونت (یعنی بدون واکسیناسیون) مصونیت ایجاد می‌شود. به 
عبارت دیگر همانطور که در شکل ۱ نشان داده شده است. محفظه عفونی تنها خروجی محفظه 
مستعد است و هیچ حالت دیگری وجود ندارد که فرد در معرض خطر به محفظه دیگری منتقل شود. 
پس از بهبودی از عفونت ناشی از بیماری» برای باقیمانده دورهه نسبت به ویروس مصون می‌شود و 
دیگر مستعد مبتلا نخواهد شد و در واقع» این یک تعریف دقیق از مدل 5115 است. از منظر نمایش 


۶۸ 


ناج 
گرافیکی همانگونه که در شکل ۱ نشان داده می‌شود هیچ ارتباطی از محفظه بهبود یافته به محفظه 
مستعد وجود ندارده یا به عبارت دیگر محفظه حذف شده حالت پایانی دینامیک بیماری است. 
این بیماری دارای دوره نهفته صفر است و این موجب می‌شود که فرد پس از قرار گرفتن در معرض 
بیماری مبتلا شود. این یک تمای زکلیدی مدل 115 از مدل 81711 است. مانند بسیاری از بیماری‌های 
عفونی» در کووید-۱۹ میانگین دوره نیفتگی بین ۴ تا ۷ گزارش شده [۱۲] که این موضوع بر پیچیدگی 
مدلسازی بمازی‌قای غقرتی مس افایی برش اه مطالنات ان داد اد که تاقلین کوند :۱۹ درراخل 
اولیه بیماری قبل از بروز علائم بالینی» موجب سرایت بیماری به دیگران می‌شوند [۶] و [۷] ۰ 
از انجایی که در مدل 115 پارامترهای 7/ و 2 ثابت در نظر گرفته شده اند» فرض می‌شود که عفونت 
زمینه‌ای در محیط های بدون مداخله خارجی تکامل می‌یابد. این یک محدودیت مهم مدل 511 
است زیرا با واقعیت متفاوت است. بسیاری از کشورها اقدامات کنترل غیر دارویی مختلفی را برای 
کاهش شیوع کوید ۱٩‏ انجام داده اند و به علاوه راه حل هایی برای اصلاح اين فرض غیر واقعی مدل 
8 در تجزیه و تحلیل داده های کوید ۱۹ پیشنهاد شده است. 
اندازه جمعیت آ به اندازه کافی بزرگ است به طوری که تعداد مواردی از جمله تعداد افراد آلوده. 
تعداد مرگ و میرها و تعداد موارد بهبود یافته را در بر دارد. به این ترتیب پارامترها در مدل 5115 را 
می‌توان به طور پایدار با دقت بالا تخمین زد. از نظر تکنیکی. این یک فرض برای مدل محسوب 
نمی‌شود بلکه شرطی برای اندازه جامعه آماری است. بنابراین از آنجایی که اين مدل در نهایت برای 
پیش بینی ریسک استفاده خواهد شد. لازم است که یک مدل مناسب با داده های قابل اعتماد تعریف 
کنیم تا علاوه بر پیش بینی . به اندازه کافی عدم قطعیت پیش بینی را نیز ارزیابی کند. 


برای درک بهتر مکانیسم بیماری‌های عفونی که توسط مدل 1 مدلسازی می‌شوند» فهرستی از ویژگی‌های 


تحلیلی این مدل ارائه می‌کنیم که دستورالعمل‌های مفیدی را برای ساختن مدل‌ها و روش‌های آماری در مدل 
۴ ارائه می‌کند. 


از آن. با اینکه سیستم دینامیکی تعریف شده توسط مدل 518 به صورت پیوسته در زمان تعریف 
می‌شود» داده های موجود در نقاط تفاتی کته و در فواصل زمانی کسسته گزارش می‌شوند. این 
بدان معناست که اگرچه مدل ۹11 دینامیک انتشار بیماری را به صورت پیوسته در زمان توصیف 
می‌کند» داده های نظارتی در مورد بیماری (مانند تعداد موارد حدید در روز با تعداد افراد بهبود یافته) 
در نقاط زمانی مشخص و در فواصل زمانی مشخص گزارش می‌شوند. این اختلاف بین مدل و داده 
۳ باید هنگام استفاده از مدل 51۳ برای تجزیه و تحلیل داده ها در نظر گرفته شود. برای مثال» اگ 
بخواهیم از مدل 51 برای پیش بینی روند شیوع نک تما ر قر ا سل استفاده کنیم باید داده ها را 


مدلسازی ریاضی و تاثیر آن در شناخت و کنترل بیماری‌ها ۶۹ 
به گونه ای اصلاح کنیم که با مدل سازگار باشند. اين کار را می‌توان با استفاده از روش هایی مانند 
فیلتر کردن یا تخمین انجام داد. بنابراین» اين یک موضوع مهم است که هنگام استفاده از مدل‌های 
ریاضی برای بیماری‌های عفونی, اختلاف بین مدل و داده را در نظر بگیریم. این کار به ما کمک 
می‌کند تا از مدل‌ها به طور دقیق تر و موّثرتری استفاده کنیم. 

و مدل 51 یک مدل قطعی است و شامل هیچ مولفه احتمالی نیست. یعنی مقادیر جمعیت در هر یک 
از گروه‌ها در هر زمان مشخص به طور دقیق توسط پارامترهای مدل و شرایط اولیه تعیین می‌شوند. 
این در حالی است که در دنیای واقعی. مقادیر جمعیت در هر گروهی با عدم قطعیت همراه است. این 
عدم فطعیت می‌تواند ناشی از عوامل مختلفی مانند خطاهای اندازگیری تغییرات در رفتار انسان 
و عوامل محیطی باشند. قابل توجه است که یک سیستم دینامیکی و یک فرآیند تصادفی دو ویژگی 
ریاضی متفاوت هستند. یک سیستم دینامیکی (به عنوان مثال مدل ۹1۳ ) لزوما تصادفی نیست. 
در حالی که یک سیستم تصادفی لزوماًپویا نیست. همانطور که در شکل ۲ نشان داده شده است» 
اندازه هر کدام از محفظه‌های (9)4, ۰ (4)[ و  17)6(‏ توابع متغیر زمان و بدون نوسانات تصادفی 
کی ور کیان میا امتهای مد وک ای ریس اه کیت ی فقو بت یزار 
زمانی که جمع آوری داده‌ها با عدم قطعیت و خطاهای تصادفی همراه است و تجزیه و تحلیل نیز بر 
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ه به راحتی می‌توان نشان داد که تعداد افراد در معرض خطر (در ورودی سیستم) یعنی (9)4, ۰ به 
طور یکنواخت غیر افزایشی است و تعداد موارد حذف شده (در خروجی سیستم) یعنی (77)0 به طور 
کشت ی کاهشی ات نوتم یشک کال مفیرخ ات انم رد هک آخاش 
که در معرض ویروس قرار گرفته اند برابر با (77)6 + (0)] < (0) 85 - ۷ بوده و به طور یکنواخت 
| 
دسته بهبود یافتگان. می‌تواند افزایشی یا کاهشی باشد. نمودار شکل ۲ به خوبی چنین خط سیری را 
که در آن زمان به اوج رسیدن و زمان کاهش به صفر (/)7 به عنوان دو نقطه عطف مهم در پاندمی 
می‌باشند. نشان می‌دهد. 

9 می‌توان نشان داد که 0 < (7)00 به این معنی است که بیماری در نهایت از بین خواهد رفت. با 
این حال. این ویژگی کاهش به صفر مشروط به مفروضاتی است که قبلا ذکر شد. نقض فرضیات 
مذکور, به احتمال زیاد باعث تداوم بیماری می‌شود زیرا یکنواختی (/)5, بیان شده در استدلال قبلی 
دیگر معتبر نیست. نمونه بارز چنین بیماری‌هایی آنفولانزای فصلی است که مصونیت در آن در طولانی 
مدت به هیچ وجه میسر نیست. 

» مدل ٩918‏ دارای ویژگی بازگشتی است. یعنی در هر زمان معین» پیشرفت بیماری فقط به مقادیر 
فعلی آنها بستگی دارد و نه به اطلاعات دیگر از گذشته. به عنوان مثال» برای پیش بینی تعداد افراد 


۷۰ م زاج 

آلوده در روز آینده» فقط به تعداد افراد آلوده در روز جاری و پارامترهای مدل 51 نیاز داریم. نیازی 
به دانستن تعداد افراد آلوده در روزهای گذشته نیست. این ویژگی بازگشتی نباید با ویژگی مارکوف 
کف متخضر | دز فرایتهای تضادفی تخت فاتون احشال شرطی استفاده شنه اسکه افتاه گر فته :شود 
خاصیت مارکوف در فرآیندهای تصادفی به این معنی است که حالت آینده یک سیستم فقط به حالت 
فعلی آن بستگی دارد و نه به حالات گذشته آن. به عنوان مثال» در یک فرآیند مارکوف سکه انداختن» 
نتیجه پرتاب بعدی سکه فقط به نتیجه پرتاب فعلی سکه بستگی دارد و نه به نتایج پرتاب های گذشته 
سکه. در مدل ,5118 هیچ قانون احتمالی در عملیات بازگشتی دخیل نیست. به اين معنی که اگر 
مقادیر اولیه جمعیت در هر یک از سه گروه و پارامترهای مدل 511 را بدانیی پیشرفت بیماری با 
قطعیت کامل قابل پیش بینی است. تمایز مفهومی بین یک سیستم دینامیکی و یک فرایند تصادفی» 
به درک تفاوت‌های بین مدل 51 و فرآیندهای مارکوف کمک می‌کند. دینامیک به معنای تغییر در 
طول زمان است» در حالی که تصادفی بودن به معنای غیرقابل پیش بینی بودن است. مدل ۹18 
یک مدل دینامیکی است. اما یک مدل تصادفی نیست. 


در طول یک بیماری همه گیر معمولاً اقدامات کنترلی مختلفی توسط دولت ها برای کاهش یا مهار گسترش 
بیماری انجام می‌شود. تأثیر مستقیم این مداخلات خارجی این است که هر دو نرخ انتقال و بازیابی» در طول 
زمان ثابت نیستند. بنابراین یک تعمیم مهم از مدل 511 » تطبیق دادن فعالیت‌های مختلف سیاست‌های 
کاهش. از جمله فاصله گذاری اجتماعی» محدود کردن حمل و نقل. پوشیدن اجباری ماسک و قرنطینه شهری 
است. همانطور که در همه گیری مداوم کوید ۱٩‏ مشاهده شد. سیاست های کاهش در طول زمان در حال تغییر 
هستند. محدود کردن جابجایی افراد مستعد و جداسازی افراد مبتلا در جمعیت. میزان ابتلا به ویروس را کاهش 
می‌دهد و منجر به کاهش نرخ انتقال بیماری یعنی (0)1/ می‌شود. در عین حال» به دست آوردن دانش بهتر در 
مورد درمان و خود مدیریتی علائم و بهبود منابع پزشکی ممکن است میزان بهبودی یعنی نرخ 7 را در طول یک 
اپیدمی افزایش دهد. بنابریان گنجاندن پارامترهای متغیر با زمان در یک مدل ۹115 منجر به گسترش مهم این 
مدل می‌شود. 

ماهیت پارامتر (0)4/ عمدتا به دو صورت قابل تعریف است. یکی این است که اجازه دهید (2)1/ یک 
تابع پارامتری باشد (مثلاًتابع نمایی نزولی) و پا یک تابع غیر پارامتری باشد [۱۴] و [۱۶]» که در هر دو 
نوع ممکن است داده های موجود تخمین زده شوند. یکی از ویژگی های مفید برای استفاده از تابع پارامتری 
7 ژمانتین است که می‌خواهیم از فصلی بودن در نرخ انتقال صحبت کنیم. در واقع این ویژگی می‌تواند برای 
مدلسازی بیماری‌های عفونی که رفتارهای فصلی دارند. مفید باشد. به عنوان مثال, بیماری‌های عفونی تنفسی 
ناشی از برخی از ویروس ها مانند کوید ۰۱۹ در برخی از ماه های زمستان با بیشترین سرعت گسترش می‌یابند. 
اين امر به دلیل عوامل مختلفی از جمله دمای پایین تر و رطوبت بیشتر در زمستان است که می‌تواند به زنده 


مدلسازی ریاضی و تاثیر آن در شناخت و کنترل بیماری‌ها ۷ 
ماندن و تکثیر ویروس کمک کند. هنگامی که فصلی بودن در نرخ انتقال گنجانده می‌شود می‌تواند به مدل‌های 
بیماری عفونی کمک کند تا رفتار واقعی بیماری را بهتر شبیه سازی کنند. اين امر می‌تواند برای پیش بینی 
شیوع بیماری و برنامه ریزی برای اقدامات کنترلی مفید باشد. به عنوان نمونه. یک تابع سینوسی می‌تواند 
برای مدلسازی تغییرات فصلی در نرخ انتقال استفاده شود. در این حالت. نرخ انتقال در طول سال با یک 
دوره مشخص نوسان می‌کند» [۱] و [۱۳] . در نظر گرفتن چنین تناوب فصلی در مدل» پیش‌بینی بلندمدت 
بهتری از یک اپیدمی ایجاد می‌کند. از آنجایی که توجه عمومی به پیش بینی همه گیری کوید ۱٩‏ به تدریج از 
کوتاه مدت به بلندمدت تغییر می‌کند. توجه به فصلی بودن اهمیت بیشتری پیدا می‌کند. در [۲] نویسنده یک 
راه ساده برای معرفی فصلی بودن ارائه است و آن اين است که فرض کنیم نرخ انتقال 0 در طول یک سال به 


صورت زير معرفی شود: 


80 - 6011 2 ات‎ ٩( (۲-16 


که در اینجا 0 متوسط نرخ برخورد و [0,1] 7 درجه فصلی بودن است که با 0  <-‏ به مدل 51۳ پایه 
تبدیل می‌شود. به علاوه (0,365] > > و در نتيجه اوج انتقال زمانی اتفاق می افتد که ع < 1. لازم به ذکر 
است که سایر توابع تناوبی يا ترکیب آنها نیز می‌توانند برای مدلسازی فصلی استفاده شوند. 


عذرمیل وا تاکسا 


فرض اندازه جمعیت ثابت در مدل ٩1‏ یک فرض ساده کننده ولی در عین حال محدود کننده است که 
ندطوز گسترده‌ای مورد انستفاده قرار می‌گیرد. این فرض به این معتی است که تعداد افراد در هر یک از سه گروه 
جمعیتی (مستعد» آلوده. و بهبود یافته) در طول زمان ثابت می‌ماند. در واقعیت. اندازه جمعیت هميشه ثابت 
نیست. افراد به طور طبیعی متولد می‌شوند و می‌ميرند. علاوه بر اين» در طول یک اپیدمی» ممکن است افراد 
به دلیل عوامل مضائی مانند مهاخرت با مرگ ناشی از بسارغه از حیعیت اف شونده. ا کته فرضی اتیازه 
جمعیت ثابت می‌تواند در درک اولیه از اپیدمی‌ها مفید باشد اما برای توصیف دقیق رفتار اپیدمی در طول زمان. 
باید دینامیک تولد و مرگ طبیعی را هم در نظر گرفت. زیرا به هر حال در یک اپیدمی طولانی‌مدت. تغیبرات 
در اندازه جمعیت می‌تواند بر روند اپیدمی تأثیر بگذارد. به عنوان مثال» در یک جامعه با نرخ تولد بالاه اپیدمی 
ممکن است سریع‌تر گسترش یابد. این به این دلیل است که تعداد افراد بیشتری برای الوده شدن وجود دارد 
ولی در یک جامعه با نرخ مرگ و مير بالاه اپیدمی ممکن است سریع‌تر فروکش کند. این به اين دلیل است که 
افراد کمتری برای آلوده شدن و گسترش بیماری وجود دارند. به علاوه اگر در یک جامعه امار مهاجرت کنندگان 
بالا باشد اپیدمی ممکن است به مناطق جدید منتقل شود. بنابراین گنجاندن دینامیک جمعیت در مدل 511 
می‌تواند به ما کمک کند تا اپیدمی‌ها را بهتر درک کنیم و پیش‌بینی‌های دقیق‌تری ارائه دهیم. با توجه به این 
توضیحات ما در ابتدا فرض می کنیم که ۸۸ نرخ تولد طبیعی و 7 نرخ مرگ طبیعی باشد. بنابراین اندازه جمعیت 
به صورت یک معادله دیفرانسیل معمولی (4) 7/۷ - (0) ۸۷ - سر تغییر خواهد کرد. در این مورد. سه 


۷۲ زاج 
خروجی برای مرگ های طبیعی وجود دارد که هر کدام در یک محفظه دح می دهد. در نتیجه مدل پایه ٩116‏ را 
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0 ۱۵ مه وا ۱6۵۵ رباص عاتا0عوگ .ظ.ل عل‌مطعصنطه اه یقع‌صهظ [1] 
,۶1141-51224 ,176 ,.1متصم‌تو؟ .لمح ععاماه 0عا عا وم ععع ی نارسمه نع 
2012 

من رعماعمط‌ن ماو مار رم مر عم ۷۵۵۱۵۵ .۰0 رتمطاببظ 2 

۰ رفت90.. م۷۷11 صطامل. :صرمل‌مصتکا 

560۵90۱۵ و ۱۱۵۵6 ۱۱6 0 قعفم ات۱6۵ 0 ۱۵/۵6۵۵ 1.1۱6 اهنا 3 

۰ ,1-15 .00 رمصتنن60 صد وامل‌م۷ آمع‌تامصمطه] ص1 

مطقصطآ۲ یا رقامتت ) رقتاماه2وواتیک ,۲.۷۷۰ بمتمصنصم؟ رت رمم‌تمن مااتاعمن رت بااعسم 4 


,1256 ,(10)6 م,.عذنا .ععلص .عم .عطلفه جم امه مام ناه که همم ۷۲۱/006 ل 
2004۰ 

,معاتتاظ ‏ ,۲.0 مجاهم ول ماه ول مت۲۵86 رآ شم۱۵ رفمتصصت رآ مجمفنم۳۵ [5] 

۰ ,4148-452 ,(442)7101 ,وان ءوع مر موه هه موه و و1۲50 

بتتقلای) پ.۲.ل ,۲۷۷۵1۲8 پل)۷۰ تالا ویک ,مق ب.ل رته۲۷۷ ویک رقتعل1 بط پلا۲۷۷ .۱0۲۲.۷ بتاضنا بکی بع۳1 [6] 
6۸0۵۲۷۱۱۰ 0 راهع عم موه وه ما ۱ 8ص 16۳7۵0۵ و کر وملا مک بصع ۷ 
0۰ ,672-675 ,(26)5 .۷۲۵0 .اع۲ 

رامآ .1.1.۲۷۷ بتتطن) باتک مصف بل.۷ معصوط یط رعصتمنا انلیا متتضا ما 1.1 بط 


آ 


۵ ۳۵6۱50۷ 0 و560 ۲۷۲۲۵ ۰ [.ظ رعصتاون ,1.۳ للملا ,..ورل مفتتو۲ .6 
مات رای ما وه ومع مس رصم اه 08۵۵0۵ ]0 
۰ ,736-742 ,(64)6 ,.عذ .101606 

0 ۵ ۱6 ۱۲۱۳۵۵۱۱۵۵0 مصحححا ,۸ ,فده انا رک بععععلط رم مصفطلط [8] 
۰ (13 ,(401 ,۲۵۲۵۲۲۲ فزن1 مکح عصهعی هه همع مان مامت 6 ۵۲ 


مدلسازی ریاضی و تاثیر آن در شناخت و کنترل بیماری‌ها ۷۳ 


نو ۶ه معط آای‌نمممطنمه مطا مه صمطتجمی ۸ .عیه یللملا , ۲۷۰۵۰ بکل‌فصنه۴ [9] 
0۰ ,700-721 ,(113)772 .90 .عص قاط اهاط وق ,عم 90 نگ ۲۵0 .تم 

2۵۲۵ ]۵ ۱۵ ۳۵۵۵006 ۰1۳6 ,تقلعل‌م؟ ظ بطاتصصه۱۱0 ۷۷ ,شش ,ماجهی .۷ پتاثناً 101 

۰ ,222021 ,(27)2 .۱۲۵0 .هگ موه متفه 10 کع همه راو وا 

بط .ات۲0 ی‌موظ فصماومی تنم وه متامصمطممط گم فصمتای‌تامرش .ع رعل‌تل‌جهاع]۷ ]11 

۹0, 44, 98-130, 95۰ 

ب0باً وبا پلالا بلاط ما ول رقطفتاا رم رقصه۷۷ بیکر ,فقا بنلط ومتای بل ,۷۷2۳8۵ نا بتاشا بش بصع 12 


۵ 0 6۱6۱۵۵9 ما ان ما ام ان 0 000ص .1 پ۲۷۳۵ .و ,۷۲۷۵ ,2 

۰ ,1915-1923 ,323 ,توف .۷۲۵ .تمحص . گ من نمرون 0017-19 
,۸۱۵۲050 ۲۰ محصامطاز ۲۷ -فمااه ریق رتطتمامطو ریق ,فمملزه)صز ۷ رظ بطمل‌همطتصاه! .۷۲۰۸۲ بنلهزهو [13] 
-(۲۷۲ 2۵ 0۲ رامع 0 مود ماه ۳60 0 رامع اه ۸۱0۵۵ 

۰ ,3550308 91 له ملطاهاته۳ .19 

و ما۵۵ وا ان 6۰۱۳0۲6 امعم رحصصون 06 ربق پهتمصتنصه [14] 
,(81)11 ,.1مظ ,اه بلاط باع۵۵0 رصع وم امه ومتعوتممصه امن عازن ۱00 
۰ ,4343-41365 

:هاگ نصا رتم۷ ها عتع یام آممتمصظ 0 ممتاصت ۲۳۵0 بط بلاط رل ,تهمنو [15] 
0۰ ,50111801 

هو 5۶۰۱۳۵ بحعط .ل باه رل بتتطر ر.۲ رعحتفتا؟ رم بجع و. تام بل مت [16] 
ب(18)3 مق هنن .ل .۵0 له اوه وه اس مان وگ 00۲۷19 [ه 
۰ ,455-482 

]0 ۲۵۵6۵ ۸ ۲۸-۰ مجمو ,۲ عهصه۷۷ رل ۲16 ربا فصقطر رو فطافموف! ۲۱ رید فمصه۲۷ ,۷ تامط2 پنآ عع1 [17] 
۰ 1 01:10) .4602-513 ,88 .تم .هو .اص عاع۱(۵0 ع عون روموت 
۰( 12402 


